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Abstract 

Quark-hadron matter at finite temperature and density is studied using a 
two-flavor, color SU(2) (P)NJL model. The hadronic effective Lagrangian, 
derived by bosonization of the quark fields and renormalized using the Eguchi 
method, emerges in the form of an extended linear a model with meson and 
diquark-baryon fields. Chiral and diquark condensates are studied as func- 
tions of temperature and baryon density. Masses of mesons and diquark- 
baryons are calculated with and without the Polyakov loop effect. We inves- 
tigate the equation of state of quark-hadron matter by taking into account 
the contributions of mesons and diquark-baryons in addition to the quark 
quasiparticles. 

Keywords: color SU(2), NJL model, lattice QCD simulation, chiral 
condensation, diquark condensation, extended linear a model, Polyakov 
loop, equation of state of quark-hadron matter 



1. Introduction 

Quarks and gluons are basic ingredients of quantum chromodynamics 
(QCD). At low temperature and density, the quarks and gluons are confined 
in colorless composites, the hadrons. However, at high temperature and/or 
high density, the quarks and gluons are de-confined and become active de- 
grees of freedom. Hence, for the discussion of strongly interacting matter in 
various situations, we ought to describe the dynamics of quarks and gluons 
and at the same time baryons and mesons in order to delineate the corre- 
sponding phases. 
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Lattice QCD (LQCD) simulations can handle strongly interacting matter 
at finite temperature and at very small baryon chemical potential. As a con- 
sequence of the sign problem in LQCD, it is extremely difficult, however, to 
deal with the QCD quark-gluon matter at finite density covering broad ranges 
of real baryon chemical potential. The Nambu-Jona-Lasinio (NJL) model 
ij is often used as an alternative, schematic approach to strongly interacting 
matter 0-0], based on chiral symmetry and its spontaneous breaking. The 
NJL model is further generalized by introducing the Polyakov loop to ac- 



count for important thermodynamical aspects of color confinement [10Hl2 



This PNJL model is by now widely used for the discussion of quark-hadron 
matter at finite temperature and density 13Nl7||. The PNJL model results 



can be compared directly with lattice QCD at finite temperature and zero 
baryon density. There is, however, a basic conceptual problem when applying 
the PNJL model at finite baryon density: color-singlet baryon formation is 
not accounted for. While color-nonsinglet degrees of freedom are suppressed 
in the "baryonic" phase by the Polyakov loop, three quarks coupled to a 
color-singlet are still delocalized and spread over all space instead of being 
confined in localized baryonic clusters. 

In this context it is instructive to study the simpler case of two-color QCD. 
With the gauge group SU(2) C , the corresponding lattice QCD approach does 
not have the sign problem, and LQCD simulations can be performed at any 
baryon chemical potential \x. Quarks carry baryon number ~. Baryons are di- 
quarks, i.e. spin singlet or triplet bosons in this theory. The resulting physics 
is qualitatively different from "real" QCD with N c = 3. Nonetheless, explor- 
ing the N c = 2 theory and designing (P)NJL type models in which both 
mesons (quark-antiquark modes) and baryons (diquarks) emerge as active 
hadronic degrees of freedom, can teach important lessons about the ther- 
modynamics and the phase structure of strongly interacting quark-hadron 
matter at non-zero baryon chemical potential. 

Several LQCD studies for color SU (2) are available [18l423| and provide 
equations of state and hadron masses of two-color quark-hadron matter at 
both finite temperature and density. It is then an interesting question to 
what extent these LQCD results can be understood and interpreted in terms 
of models, such as the PNJL approach, in order to identify leading mech- 
anisms and basic symmetry breaking patterns. Our work is aimed in that 
direction. It presents major steps beyond a previous NJL model calcula- 
tion 2m performed in a similar context, and it complements other related 



two-color QCD studies [25H39 
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In this paper, we study the color SU(2) NJL model using bosonization 
and renormalization techniques to construct the meson and diquark-baryon 
effective Lagrangian. The derivation of the complete hadron Lagrangian 



follows the Eguchi method |40j and renormalization lead to a generalized 
linear a model with Pauli-Giirsey symmetry. The partition function and 
the thermodynamical potential are derived with and without incorporation 
of Polyakov loop effects. The equation of state is deduced and the phase 
diagram is studied in detail. We investigate the excitation spectrum of quark- 
hadron matter at various temperatures and densities. We then include the 
effect of hadrons in the thermodynamical potential and calculate the pressure 



and the quark density using the Gaussian approximation (4l| of the hadron 
Lagrangian, taking only the mass terms into account. 

This paper is organized as follows. In Sect. 2, the NJL model for the 
case of color SU (2) with two flavors is introduced, bosonization is performed 
and an extended linear a model is derived. The thermodynamical potential 
of the meson and diquark-baryon fields is deduced in Sect. 3. In Sect. 4, we 
present numerical results and compare with LQCD results. Sect. 5 is devoted 
to conclusions and outlook. 

2. Two color NJL model 

2.1. NJL Lagrangian and auxiliary fields 

We start with the NJL Lagrangian for color SU(2) and two quark flavors: 

G 

£njl = 4>{ilvd u - m + 70/i)^ + -^-[{ipip) 2 + {ipijsrip) 2 ) (1) 
+ ^(^t2T 2 Cij T )(ij T Ct l5 t 2 T 2 ^) . 

The quark fields are ip = (u, d) T , with masses m u = = mo- The trans- 
pose of ip is denoted by tp T . The SU(2) isospin and color matrices are 
written n and tj, respectively. We introduce the quark chemical potential 
\x for the study of finite baryon density, with the baryon chemical potential 
fis — 2{i. The diquark-baryon fields involve the charge conjugation operator 
for fermions, C = 27072, together with t 2 and t 2 to form color and isospin 
singlets. Starting from a color current-current interaction proportional to 
(ip r Yi 1 tjip)(ijj'yL 1 tjip), a Fierz transformation provides equal coupling constants 
Gq for the meson channels and H for the diquark channels, featuring the 
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Pauli-Giirsey (PG) symmetry. We keep, however, the coupling constants in- 
dependent so that we are able to study as well cases in which PG symmetry 
is not exactly realized. 

The bosonization technique is now used to write the Lagrangian in terms 



of auxiliary meson fields, a(x), 7?(x), and diquark fields, A(x), A*(x) |40l. 142 



44j . The partition function of the NJL model including these auxiliary fields 



is 

Z = J V^V^VaVnVAVA* exp (i J d A xC au ^j , (2) 

with the Lagrangian 

C - ib T C)( g d i 5 t 2 T 2 A{x) \ f if> Von 

Caux ~ 2 ( ^ V C >{ -g d A*(x) l5 t 2 r 2 ^(-/i) ) { C$ T )&> 

~\M 2 Aa\x) + n\x)) - 1 -M 2 d A*(x)A(x) , 

with S' _1 (±/i) = i^^d^— m ±7o/i — go(a(x)±i^ 5 7i(x)-T). We have introduced 
a meson coupling constant g and a diquark coupling constant together 
with bare scalar and diquark masses, M s and M d , in preparation of the 
standard renormalization scheme. The meson and diquark coupling constants 
and masses are related as jfa = G and jfr = H . 

Introducing a chiral order parameter a"o (proportional to the quark con- 
densate (^V)) as a { x ) = °~o + s(x) and the diquark condensates, A and Aq, 
as A(x) = A + d(x) and A*(x) = Aq + d*(x), and integrating out the quark 
fields, the effective Lagrangian entering the partition function (j2j) becomes 

C eff = --tr(ln£- 1 + ln(l + £A')) (4) 
~ \mI{s\x) + 7f 2 (x)) - M*a s(x) 
-^MjA*A - l -Ml{Ald{x) + d*(x)A ) - ±M 2 d d*(x)d(x) . 

The trace is taken over spin, flavor and color spaces. The matrices S~ l and 
K are defined as, 

o-i f S \n) A \ 



ft _ ( -9o{s{x) +ij 5 T ■ 7f(x)) 9dJshr2d(x) \ , g x 

V V -9dd*(x)-f 5 t 2 r 2 -9o(s(x) -ij 5 r • 7r(x)) J ' 

with A~ = g d j5hr 2 A , A + = -gd^llshri and Sq 1 {±h) = «7 M 9 M - m± 70/x. 
A dynamical quark mass is denned as m = mo + go&o- The Nambu-Gorkov 
quark propagator matrix S is determined by solving SS~ X = 1 and expressed 
as: 



with the components 



G± = (S^i/O-E*)" 1 (8) 
E± = A^5 (T^)A ± 
^ ± = -5 , o( T/U )A ± G ± . 

A simple form for the components of the Nambu-Gorkov propagator is 
found introducing the energy projectors onto states of positive and negative 



energy for free massive spin 1/2 quasi-particles following Huang et al. 45l. 146 



with = yp 2 + m 2 . These operators satisfy the projection properties 
A±A± = A±, A±A T = 0, A + + A_ = 1, and similar relations for A. Fur- 
thermore, the two projection operators are related as 7cA±7o = A T and 
7 5 A±7 5 = A±. The Dirac propagator can be written 

Po - •C'f PO + -frp 

Taking its inverse we arrive at 

^(ip) = (po - ^J)7oA+ + (po + ^)7oA_ , (11) 

with E^ 1 = E p ± ft. It follows that the components of the Nambu-Gorkov 
propagator are written as 

Po - (#a) - (^a) 

A ± A ± 

H = ~2 7^±Y2 A + + ~2 / FT N2 A - ' ( 13 ) 

Po - (^a) 2 Po - (^a) 2 
5 



The quasi-particle energy is = * /(E±) 2 + g\\ A | 2 with the dynamical 



quark mass m = m + g a and the diquark gap A . 

2.2. Mean field approximation 

Keeping only the expectation values (cr(x)) = oq, (A(x)) = A and 
(A* (a;)) = Aq, i.e. dropping the fluctuating meson and diquark fields s(x), 
tt(x), d(x) and d*{x), one arrives at the mean- field Lagrangian 

Cmf = -^trlnS- 1 - l -M 2 s al - hl 2 d A* A . (14) 

The first term is calculated using the relation trlnS 1-1 = lndetS*" 1 . The 
determinant is taken over spin, flavor, color and momentum space and calcu- 
lated using Eq. ( ITT]) together with the properties of the projection operators 



45, 46 . One finds 



det S- 1 = ^/ (p% - (E^)( P l - (E+)i) . (15) 
Temperature is now introduced using the Matsubara formalism, 

with the replacement P o — > iu n where u n = (2n + 1)ttT are the fermionic 
Matsubara frequencies. Taking the frequency sum, the thermodynamical 
potential becomes 



-2tr J ^T[ln(l + e-^) + ln(l + e-^A)] 



+\m 2 A + \mI\a^ . 

Here we have written AqAo = |Ao| 2 and the inverse of the temperature T as 
f3 = 1/T. 

The derivatives of the thermodynamical potential (1171) with respect to a 
and |A | determine the chiral condensate and the diquark condensate at the 
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minimum of £Imf : 
OQmf 



da 



-tr 



d 3 p gom 
(2vr)3 E p 



Et 



E: 



E 



^ (1 - 2f(E+)) + ^ (1 - 2f{EZ)) 



+M> = 



lAoltr 



d 3 p 
(2tt) 3 



-L (l - 2/(25+)) + -L (1 - 2f{ED) 
E A E A 



+M 2 d \A \ = 



(19) 



with the Fermi distribution function f(E) = (1 + e l3E ) 1 . The chemical 
potential /i is included in the definition of E A . 

2. 3. Hadron Lagrangian from bosonization 

We would like to extract the properties of mesons and diquark-baryons at 
finite temperature and density. For this purpose, we expand the logarithmic 
term in Eq. (j3J) as 



- -tr\n(l + SK) 
2 



k=l ' k=l 



(20) 



where 



i(-l) 



fc+i 



2k 



-tr{SK) k . 



The matrix 



SK 



A B 
C D 



has elements 

.4 
B 
C 
D 



-G + g (s(x) + ij 5 T ■ 7? (x)) - H~g d d*(x)^f 5 t 2 r 2 
G + g d r y5t 2 T 2 d(x) - H~g (s(x) - ij 5 f ■ tt(x)) , 
-H + g (s(x) + ij 5 T ■ n(x)) - G~g d d*(x)^ 5 t 2 r 2 
H + gdl5t 2 r 2 d(x) - G~g Q (s(x) - i^ 5 r ■ vf(x)) . 



(21) 



(22) 



(23) 



With these expressions we are able to derive the Lagrangian involving the 
meson and diquark-baryon fields. 
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-trSK , to derive the 



2-4- The gap equations (k=l) 

We first work out the case of k = 1, that is U^ l > 
mass gap equations, which correspond to the minimization condition f|T8|) of 
the thermodynamical potential in the mean-field approximation described in 
the previous subsection. 



2 tr 



^TTU [do(-G + -G )s(x)+ g (-G + i>y 5 + G r^f-frl 



(2vr) 4 

-g d H~j 5 t 2 T 2 d*(x) + g d H + j 5 t 2 T 2 d(x)] 

/^(r, sW + r idW + r, d . (l 



(24) 



The trace of the Dirac matrix in the pion term gives zero and we have dropped 
this term here. The result for I\ is 



:9o tr 



d A p 



(G + + G ) = 2img tr fc 



d 4 p 



x 



(25) 



IpI - (Elf pi - (E A y E p \pl - {Elf 

The trace tr is taken over spin, flavor and color spaces and try c is over flavor 
and color only. This expression agrees with the mean-field equation (fT8l) at 
T = 0, when explicitly writing = E p ± /i, so that 



on 



MF 



dan 



-9o m tr 



The momentum space integrals (1251) are divergent and must be treated 
accordingly. For the second term under the integral we write: 



d 3 p 



-L (1 - 2f(E+)) + -L (1 - 2/(££)) (26) 

.-^A -^A 



1 



+ M s 2 a = . 



JL 

E p 



vl 



(El) 2 



vl 



(EX 



(27) 



e p \vl - (E^y 



i' 



E P [p 2 - (EX 



2a z 



■4:uE p 
1 + 



Pa 



(E+y + 4iiE, 



AuE p 



vl 



(Pi - (El 



+ \2\2 



+ 



(£a) 2 

1 



p\ 



(Ely 



AuE p 



pI 



(Ely 



(pi - (EX) 2 ) 2 
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where in the third line only those terms have been kept that are divergent 
when taking the momentum integral. Using this expression T s can be written 

as 

T s = 2mg {I 2 - 2/i 2 / ) , (28) 
with the divergent integrals I 2 and Iq defined as follows: 




The result for T d is 

T d = ~g d tr J ^H+ l5 t 2 r 2 = g 2 d A* I 2 , (31) 

while 

IV =g 2 d A I 2 . (32) 

At this stage the effective Lagrangian is C e ff = Cmf + £^ with 

= 2g m(I 2 - 2/2 2 I )s + g 2 d I 2 (A* d + A d*) . (33) 

It involves the leading-order couplings of the condensates to the meson and 
diquark fluctuations, s(x), d(x) and d*(x). 

T s together with other terms linear in s in the effective Lagrangian gives 
the mass gap equation by requiring 

- M> + 2g m(I 2 - 2fi 2 I ) = . (34) 

We recall that the relation between m and 00 is m = g <r + m . With the 
T d and IV terms and the linear diquark-baryon terms in the Lagrangian, we 
obtain the diquark condensate equation, 

- M 2 A + 2g 2 A I 2 = . (35) 
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For later use, we write I2 and Iq as functions of temperature and chemical 
potential using the Matsubara formalism: 



tr/. 



n 

d 3 p 



d 3 p 



+ 



1 



{ED 



(36) 



(2vr)3 \2El 



1 .(l-2/(£+)) + ^(l-2/K)) 



2£; 



Proceeding analogously in Eq. fl25|) . the mass gap equation is modified 
accordingly and agrees with Eq. fl26|) . 
Next, we work out Iq\ 



tr 



tr 



<i 3 p 

(2^p V((iw„) a -(^) 2 ) a ' ((^n) 2 - (£ A ) 2 ) 2 
4^53 ( X - 2 ^a) - 2£+/3e^/ 2 (£+)) 



n 

d 3 p 

l— 3 (l-2f(E A )-2E^(3e^f(E A ] 



(37) 



4(^a 

With these expressions, the gap equations are - as expected - identical to 
the expressions found by taking derivatives of Omf with respect to the con- 
densates (see Eqs. ffT8l) and ffT9|) ). The divergent integrals Jo and I2 are 
calculated introducing a momentum cut-off A. 

2.5. Mass and kinetic energy terms (k=2) 

The kinetic energies and the mass terms of mesons and diquark-baryons 
emerge from Eq. fT2~TT) at the order of k = 2: 



(38) 



= - A ti{skf = I d 4 xd 4 y J2^-y)^)My) > 



where 0« = a, 7?, d, d*. The non-local propagator terms are expanded as 



(39) 



S(x - y)T(0) - id»5(x - y)d fl T(q) \ q=0 



--d^d v 8(x-y)d l ,d u T(q) 



9=0 
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Detailed results for the quantity T(q) \ q=0 and the derivative terms d fl T(q) \ q=0 
and d ll d v Y{q) | g= o are given in Appendix A. The emerging piece of the effec- 
tive Lagrangian at order k = 2 is: 

£ (2) = l -glh{d^f + g l{I 2 -2I^ (40) 
+ lg 2 I (d,s) 2 + g 2 (I 2 -2I (^ + m 2 ))s 2 
+ l -g%(d»d*)(d,d) + g 2 (I 2 - I g 2 \A \ 2 )d*d 

- tg 2 I fi(d*d d - dd d*) - \g%{A* 2 d 2 + A 2 d* 2 ) 

- 2g g 2 I ms(A* d + A d*) . 

The Lagrangian density has mass dimension 4, and Jo an d I2 have mass 
dimensions and 2, respectively. Hence when I appears in the coefficients 
there are factors of mass dimension 2 such as squares of the derivative <9 M , the 
chemical potential /x, the quark mass m and the gap energy Ao and Aj$. The 
pion term has a simple form due to the pseudo-scalar nature of the pion field, 
while the scalar boson s(x) couples to the diquark-baryon fields through the 
diquark condensates. 

2.6. The coupling terms (k—3) 

The coupling vertices involving combinations of three mesons or diquark 
fields are obtained from the third-order term: 

C/( 3 ) = -Ltr(SK) 3 = -^tr{A 3 + ABC + BCA + CAB (41) 

+CBD + BDC + DCB + D 3 ) . 

Non-local terms with derivatives of hadron fields vanish, and is con- 
structed with only local terms. The details of the derivation are found in 
Appendix A. The final result for the k — 3 part of the effective Lagrangian 
is: 

£ (3) = -2g I m[g 2 (s 2 + 7r 2 )s + g 2 d d*ds] (42) 
-9po His 2 + n 2 )(A d* + A* d) + g 2 d*d(A d* + A* d)} . 

The coupling terms have mass dimension 3 due to the presence of altogether 
three meson and diquark-baryon fields. Hence, the coefficients of all k — 3 
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vertices have mass dimension one with factors such as to, A and Aq, together 
with the dimensionless integral Iq. Those terms with to in the coefficient 
conserve the baryon number, while terms with Ao or Aq are baryon number 
changing. 

2.7. The interaction terms (k—4) 

Four-point interaction terms involving mesons and diquarks derive from 
the fourth-order term, 

f/W = ^tr(S7i~) 4 . (43) 
8 

The non-local pieces vanish. Details of the derivation are again given in Ap- 
pendix A. The final expression for the k = 4 part of the effective Lagrangian 
is: 

= -g 4 I s 2 n 2 -glg 2 d I (s" + n 2 )d*d (44) 
-\gth{s' + ^)-\ g %d^£ 

= -\h[gl^ + ^) + g ld*d] 2 . 

It generates interactions between diquarks, scalar mesons and pions in all 
possible combinations. 

2. 8. Identification with a generalized linear a model 

The limiting case of exact Pauli-Gursey symmetry in color 577(2) is re- 
alized with g — 9d and M s = M d . The fourth order Lagrangian density is 
written as 

C^ = - l -gMs 2 + K 2 + d*dy . (45) 

This form and the structure of the other terms suggest, not surprisingly, 
that the hadron Lagrangian is related to a generalized linear a model. To 
demonstrate this, we introduce a fourth-order extended a model Lagrangian 
density with a, 7?, A and A* fields: 

4 4) = -^oW + ^ 2 + A*A) 2 . (46) 
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Separating again mean fields and fluctuations by the relations a = cr + s(x) 
and A = A + d(x) together with its complex conjugate, one finds 

a 2 + vf 2 + A*A = (a + s(x)) 2 + n 2 (x) + (A + d(x))(A* + d*(x))(47) 
= a 2 + A A* + 2a s(x) + A d*(x)+A*d(x) 
+s 2 (x) + 7? 2 (x) + d*(x)d(x) , 

and the 4-th order Lagrangian (|46j) is rewritten as: 

4 4) = ~g^oH + 2v 2 (2a o s + A d* + A*d) (48) 

+2v 2 (s 2 + tt 2 + d*d) + Aols 2 + 2\ A \ 2 d*d 
+A 2 d* 2 + A* 2 d 2 + Aa A sd* + Aa A*sd 
+4(T (s 3 + sd*d + svr 2 ) + 2A (d* 2 d + s 2 d* + vr 2 cf ) 
+2A* (d*d 2 + s 2 d + 7? 2 d) + (s 2 + vf 2 + d*d) 2 ] , 

where we have defined a 2 + |A | 2 = v 2 . 

Hence one can write the hadron Lagrangian in the compact form of a 
generalized linear a model with inclusion of diquark fields: 

Ce„ = \glh[(d^) 2 + {d u a) 2 + \(d u -2in5 uQ )A\ 2 ] (49) 
- 1 -{M 2 - 2g 2 I 2 + Ag 2 I 0f x 2 )(a 2 + tt 2 + A*A) 
-i, %[(a 2 + vf 2 + A*A)-, 2 ] 2 . 

In order to write the Lagrangian in this form, we have added —2g 2 I fi 2 A*A 
to the mass term and combine with the a 2 + tt 2 terms to arrive at the compact 
symmetric mass term shown in the second line. We then add the counter 
part, +2g 2 : I fi 2 A*A, to the kinetic energy term proportional to (d u A*)(d u A) 
and arrive at the compact form ~g 2 I \(d v — 2ifi5 u0 ) A \ 2 . 

Consider next the explicit chiral symmetry breaking term. Obviously, we 
have a term linear in the scalar field, 2(yf ^o(-^2 — 2/z 2 / )<s, which is proportional 
to the bare quark mass. This suggests the explicit chiral symmetry breaking 
term of the form 

C S b = 2g m (I 2 - 2/i 2 / )a . (50) 
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This term provides a new mean-field mass equation including explicit chiral 
symmetry breaking, 



- M S V + 2(7o(m + g a ){I 2 - 2/j 2 I ) = . (51) 

Self-consistent solution of this gap equation determines the dynamical quark 
mass, m = m + #oCo, that appears also in the integrals I2 and Iq. 

2. 9. Renormalization 

The hadron Lagrangian ( l4"9"j) extracted from the NJL model involves the 
divergent quark loop integrals I2 and Jo- The NJL approach is valid for quark 
momenta below a characteristic scale, A ~ 0.6 GeV, at which the integrals 
are cut off in practice. Nonetheless, when written in the form (l4"9|) as a 
generalized linear sigma model, renormalization has to be performed. We 
follow the Eguchi method [40] . The kinetic term implies the wave function 
renormalization, 

glh = zil (52) 

1 1 
a = Z M a R 7? = Z^TTr 

A = ZIAr A* = zIa* r . 

The mass renormalization condition is 

M 2 S - 2glh + Aglhn 2 = . (53) 

The coupling constant renormalization is written as 

2^ 4 / = Z-'Xo, (54) 
A = Z 2 M Z X 1 \ (s . 

The mean field is renormalized as 

1 

v = ZIjVq , (55) 
and the explicit chiral symmetry breaking term is renormalized as 

2g m (I 2 -2fi 2 I ) = Z s l B e (56) 
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When expressed in terms of renormalized fields and couplings, the effec- 
tive Lagrangian (I4"9"|) reads: 

C e «sb = ^[(d u n R ) 2 + (d u a R ) 2 +\(d u -2tfi5 u0 )A R \ 2 } (57) 

-\{°r + 4 + A fi A « - ^ 2 ) 2 + ^ . 

Together with the quark term, we have now a consistent Lagrangian to de- 
rive the thermodynamical potential for the interacting quark-hadron system. 
The hadron dynamics itself is governed by the generalized linear a model La- 
grangian fl57|) . The explicit chiral symmetry breaking term, ea R , is related 
to the bare quark mass mo- When the PG symmetry is not satisfied, the 
meson part and the diquark-baryon part are renormalized independently. 



3. Thermodynamics of mesons and diquark-baryons 
in the Gaussian approximation 

We now discuss the contributions of mesons and diquark-baryon fields 
to the thermodynamical potential, using the Gaussian approximation. The 
starting point is the extended linear a model Lagrangian derived in the pre- 
vious section. The whole Lagrangian is written as 

£ = C>MF + £hadron ■ (58) 

The mean field part has been worked out previously to provide the thermody- 
namical potential Omf- The extended linear a model includes hadron fields 
up to fourth order. The complete integrations over the hadron fields cannot 
be performed to this order. In the present work we restrict ourselves to the 



Gaussian approximation 41] , taking only the second order (mass) terms into 



account and perform the hadron integrals for the partition function as 

Z = Z q J VsVnVdVd* exp J d A xC {2) adro )j , (59) 
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with 

Z q = exp ^-^Qmf^ , (60) 

+-(d"d*)(d tl d) - ^m 2 d d*d - ifi{d*d d - dd d*) 

--A* 2 d 2 - -A 2 d* 2 - 2mA* sd - 2mAsd* . 
2 2 

Here we have used the renormalized fields but omitted the index R for 
simplicity. For further convenience, we have introduced A = and 
A* = (^Aq in the above equation. The renormalized masses are defined 
as 

2 



in 

,2 r 7\ /r2 



m s 
m 2 



[M 2 - 2g 2 (I 2 - 2/ 0/ , 2 )] /(g 2 I ) , (61) 
[M 2 - 2g 2 (I 2 - 2/ (/i 2 + m 2 )] /(g 2 I ) , 
[M 2 -2^(/ 2 -/ |A| 2 )] /(g 2 I ). 

The meson and diquark-baryon masses flBTl) found by bosonization and renor- 
malization are identical to the masses obtained by solving the corresponding 
Bethe-Salpeter equations taking only the divergent integrals Iq and I 2 in the 
quark-loop integrals. Of course, using the Gaussian approximation means 
ignoring important interaction terms generated by C^ 3 ' and C^ 4 \ These will 
be dealt with in separate forthcoming work. A discussion of some qualitative 
features of such terms will follow the concluding section. 
The hadronic partition function is 

Zhadron = j VQVtT exp (-S®) , (62) 

S {2) = l ~Y. f d 3 x[<5>\x)M n (Q)<S>(x)+n\x)N n (Q)n(x)} , 
Q J 

where Q = (icu n ,q), with ui n = 2irnT being the boson Matsubara frequency 
and J2q = TJ2nJ (f^p ^ ne fi^ds ^ are defined as $ T = (s,d, d*) and 
tt t = (7r + , 7T°, The matrix M n (Q) is the mass matrix in Euclidean space 
with Q 2 = q 2 + u 2 and Qo = ioj n : 

/ Q 2 + m 2 2mA* 2mA \ 

M n (Q)=\ 2mA l(Q2 + m 2)_ 2/i g o A 2 (63) 

2mA* A* 2 UQ 2 + m 2 ) + 2fiQ l 
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On the other hand, the matrix N n (Q) of the pion term is diagonal and 
proportional to the identity matrix as 



The hadronic thermodynamical potential, Qhadron 
comes 



(64) 

"t7 ^ ^ hadroni be- 



n 



hadron 



~ ( lndet M n{Q) + Indet N n (Q)) 



(65) 



+ 3( y + Tln(l-e-' 3 ^; 



^/q 2 + mf. The excitation energies 



The pion excitation energy is cu^g) 
o;j(i = 1,2,3) are the eigenvalues of the matrix M(u,q) in Appendix B 
obtained by solving the dispersion relation det M(co, q) = with respect to 
u numerically. The momentum integrals of these boson zero point energies 
are divergent. We regularize the integrals by the cut-off momentum A as in 
the case of the quarks and subtract the zero point energies in the vacuum, 
T = fi = 0. We can verify that the zero point energy contribution van- 
ishes when |A| =0. At finite |A| there appears a finite contribution to the 
thermodynamical potential. 

A related study of the thermodynamical potential beyond the mean field 
approximation, including hadron contributions in the two-color NJL model 



by using the Gaussian approximation, has been performed by He |37J. The 
effect of the pair mode in condensed matter has been investigated by Diener 
et al. 
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This method is applied to the NJL model and the effect of 
the hadrons is estimated by a perturbation method for the quark density 
An expansion of the N c = 3 two-flavor PNJL model beyond mean field 
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approximation has been reported in ref. [14 . 



4. Numerical results 

4-1. Polyakov potential 

For the discussion of the equation of state of quark-hadron matter, it is im- 
portant to take into account the quark confinement effect. We adopt here the 
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Fukushima method and add the Polyakov loop potential [/($ [A], <&*[A]; T) 



to the two-color N JL Lagrangian . The derivative <9 M acting on the quark 
field is replaced by the covariant derivative = — iA^. The temporal 
background color gauge field is introduced A 4 = iA° with A = gA®^ and 
the SU(2) Pauli matrices t a (a = 1,2, 3) in color space. In the Polyakov gauge 
this temporal gauge field is diagonal in the color space. For the SU(2) color 
group it is represented as t 3 9 where 6 is real. 



The Polyakov loop potential U is written as [32 

[/($, T) = -6T(24$ 2 e^ a + ln(l - $ 2 )) , (66) 
in terms of the traced Polyakov loop 

$ = ^-tre i/M4 = cos(/3fl) . (67) 

The log term comes from the Jacobian in the 577(2) color space. The two 
parameters a =858.1MeV and 6 1/3 =210.5MeV are taken from the work of 
Brauner et al. j32|. 

Taking into account the Polyakov loop effect implies the replacement 

ln(l + e-^) ^ I i n (i + 2$ e "^ + e" 2 ^) . ( 68 ) 

Hence, the Fermi function f(E) in the mean field equation is replaced by 

1 + §e l3E 

m = l+2 ^ E+eWE ■ (69) 

This means that single quark quasiparticles are suppressed in the hadronic 
phase as $ 0, whereas color-singlet combinations with two quarks survive. 

We now proceed to discuss numerical results. For this purpose we fix the 
parameters of the PNJL model following the paper of Brauner et al, the 



list of which is given in Table I of 32|. These parameters are set such that 



the pion mass in free space is m, = 140 MeV and the pion decay constant 
is /„ = ^/2/3 x 93MeV=76MeV, reflecting the ^/W c behavior of f n . The 
coupling constants are G = H = 7.23 GeV -2 , the cut-off momentum is 
A = 657 MeV and the bare quark mass is m = 5.4 MeV. 
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4-2. Phase diagram and order parameters 

We start by showing the chiral condensate, the diquark condensate and 
the Polyakov loop with and without the Polyakov loop coupling effect. These 
condensates are obtained by solving the gap equations for the chiral conden- 
sate <7 and the diquark condensate |A |. Results are presented in Figs. 
[T] and [2] as functions of temperature at zero chemical potential (/i = 0). 
The crossover temperature of the chiral condensate defined at a c = |cr 




100 200 300 400 50' 



T [MeVl 




100 200 300 400 50' 



T [MeVl 



Figure 1: The chiral condensate oq in 
unit of ctq(T = 0) at zero tempera- 
ture (solid curve) and the Polyakov 
loop <I> (dashed curve) as functions 
of temperature at zero chemical po- 
tential (/U=0) for the case of the NJL 
model with decoupled Polyakov loop. 



Figure 2: The chiral condensate o"o 
(solid curve) and the Polyakov loop 
<I> (dashed curve) as functions of tem- 
perature at zero chemical potential 
(/z=0) for the case of the PNJL model 
including coupling of quarks to the 
Polyakov loop. 



is 180 MeV. The critical temperature of the Polyakov loop is set at T = 
270 MeV for the pure gauge case without coupling to quarks. At this point 
we adopt the N c = 3 value of To as in ref. 32j. According to ref. [42], the 
^-dependence of T goes approximately as T /i/a ~ 0.6 + 0.45/iV^, sug- 
gesting roughly a 10% difference between T (N C = 2) and T (N C = 3) that 
we can ignore. When the coupling between the quark and the Polyakov loop 
is introduced in the PNJL model, the chiral and deconfinement crossover 
transition temperatures become strongly correlated as shown in Fig. The 
crossover temperatures of both order parameters defined at half of their full 



values are now T v ~ Ta 



225 MeV. 



- x — dec 

Results at finite chemical potential (fi = 100 MeV) are shown in Figs. 
E] and HJ In this case the diquarks condense and the chiral condensate is 
reduced smaller where the diquark condensate is finite. Again the coupling 
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Figure 3: The chiral condensate do 
(solid curve) , the diquark condensate 
|Ao| (dash-dotted curve) in unit of 
ao(T = 0) and the Polyakov loop 
(dashed curve) as functions of tem- 
perature at finite chemical potential 
(/x = lOOMeV) for the case of the 
NJL with decoupled Polyakov loop. 



Figure 4: The chiral condensate o"o 
(solid curve) , the diquark condensate 
|Ao| (dash-dotted curve) in unit of 
(Jo(T = 0) and the Polyakov loop $ 
(dashed curve) as functions of tem- 
perature at finite chemical potential 
(H = lOOMeV) for the case of the 
PNJL model including coupling of 
quarks to the Polyakov loop. 



of quarks to the Polyakov loop makes the deconfinement transition smooth, 
so that both chiral and deconfinement crossover temperatures become quite 
similar. 

The phase boundary of the de-confinement transition ($ ~ 0.5) is insen- 



sitive to the chemical potential as discussed by Brauner et al. [32(. Hence, 
we plot in Fig. [5] and [6] only o"o and | Ao | for various chemical potentials as 
functions of temperature. The chiral condensate <To shown by the top smooth 
curve (/z = 0) in Fig. \5\ stay unchanged until diquark condensate sets in at 
fi = \ra^ = 70 MeV. ctq at \i = 75 MeV is depleted in the small temperature 
region, where the diquark condensate |Ao| is finite as shown in the right hand 
figure. This behavior continues as the chemical potential increases as shown 



for ctq and |Aq| at fj, — 100 and 200 MeV |24j. |26I|. These behaviors agree with 



the results shown in Fig. 3 of ref. [32 



4.3. BCS-BEC crossover 

The two gap equations ( !T8|) and (TT9T) for a Q and | A | point to an interesting 
interplay concerning the BEC-BCS crossover. When the chemical potential fi 
exceeds the critical point of diquark condensation, the dynamical quark mass 
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Figure 5: The chiral condensate o"o in 
units of a (T = 0) for /i = 0, 75, 100 
and 200 MeV (from top) as functions 
of temperature. 
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Figure 6: The diquark condensate 
|Ao| in units of (Jq{T = 0) for fj, = 0, 
75, 100 and 200 MeV (from bottom) 
as functions of temperature. 



m = go<jQ + m is much larger than the chemical potential fi as shown in Fig. 
[7J Hence, the system is in the BEC phase. This is demonstrated in Fig. [8] 
showing the quasi-particle energy of quarks, E^(p) = ((E~) 2 + gj| AqI 2 ) 1 / 2 , 

with E~ = \/ m? + p 2 — fi at fi = 100 MeV. At this chemical potential the 
minimum of the quasi-particle energy is at zero momentum. As fi increases, 
the gap | Ao| increases, while the dynamical quark mass m decreases as shown 
in Fig. [7J Once the chemical potential /i exceeds the dynamical quark mass, 
H > m, the system undergoes the BEC-BCS crossover and turns into the BCS 
phase as /x increases further. The quasi-particle energy E^ at ji = 200 MeV 
shown in Fig. [S] has a minimum at finite momentum. The crossover point 
is insensitive to the temperature (see Figs. and E]) as long as the diquark 
condensation is finite. 

4-4- Meson and diquark masses 

We start with scalar sigma boson and pion masses at zero chemical po- 
tential. In this case, the diquark condensate is zero and there is no sigma- 
diquark mixing. The behavior of these masses is similar to that found in the 
SU(3) color NJL model. The only difference is the role of the Polyakov 
loop for two colors here as compared to the standard color SU(3) case. 
We show the sigma and pion masses as functions of temperature with and 
without the quark-Polyakov-loop coupling at zero chemical potential (/i=0) 
in Figs. [9] and [KB The explicit chiral symmetry breaking by the small 
quark mass, ttlq = 5.4 MeV, gives the Nambu-Goldstone pion a small mass, 
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Figure 7: The dynamical quark mass 
m (solid curve) and the chemical po- 
tential \x (dashed curve) as functions 
of the chemical potential. The BEC- 
BCS crossover point corresponds to 
the crossing point, \x = 120 MeV. 
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Figure 8: The quasi-particle energy of 
quarks at \i = 100 MeV (solid curve) 
and at (j, = 200 MeV (dashed curve) 
as functions of momentum p. 



m n = 140 MeV at zero temperature, while the sigma mass, m s = 610 MeV 
stays around twice the dynamical quark mass. For zero chemical poten- 
tial, the diquark-baryon masses are degenerate with the pion mass due to 
the PG symmetry. As the temperature increases, the pion mass starts to 
increase, while the sigma mass decreases until the temperature approaches 
the crossover temperature of about 180 MeV. The two masses meet at the 
crossover temperature and increase jointly as shown in Fig. [9j 

When the quark-Polyakov-loop coupling is introduced, the pion mass in- 
creases earlier with temperature as shown in Fig. HOI The reason is that the 
quark mass drops more slowly with increasing temperature as shown in Fig. 
|2j The sigma mass stays almost constant up to the crossover temperature 
and then both the pion and sigma masses increase rapidly with increasing 
temperature. 

In the case of diquark condensation at finite chemical potential /i > 
all scalar hadronic degrees of freedom have non-diagonal coupling terms. The 
mass matrix therefore needs to be diagonalized together with the chemical 
potential terms. We discuss the behavior of the masses after diagonalization 
in Appendix B. The renormalized pion mass is written as 

ml ={M^-2gll 2 + Agll^)l{gll ) = f Q . (70) 

The pion mass stays constant at 140 MeV and proportional to the square 
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Figure 9: Sigma (dashed curve) and 
pion masses (solid curve) in MeV as 
functions of temperature in MeV at 
zero chemical potential (/U=0) for the 
case of the NJL model with decou- 
pled Polyakov loop. 



Figure 10: Sigma (dashed curve) and 
pion masses (solid curve) in MeV 
as functions of temperature at zero 
chemical potential (/U=0) for the case 
of the PNJL model with coupling of 
quarks to the Polyakov loop. 




Figure 11: The sigma (dotted curve), 
the pion masses (solid curve) and 
diquark masses (dash-dotted and 
dashed curves) in MeV as functions 
of the chemical potential /i at zero 
temperature (T = 0) . The sigma and 
diquark coupling is dropped when 
|A |/0. 



Figure 12: The sigma (dotted curve), 
the pion masses (solid curve) and 
diquark masses (dash-dotted curve) 
and (dashed curve) in MeV as func- 
tions of the chemical potential /j, at 
zero temperature (T = 0). The 
sigma and diquark coupling is intro- 
duced when |Aq| / 0. 
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root of the bare quark mass until the onset of diquark condensation. The 
sigma mass is larger than the pion mass, m 2 = m 2 + 4m 2 , with m the 
dynamical quark mass as discussed in Appendix B. Meson and the diquark- 
baryon masses are shown as functions of the chemical potential fi at zero 
temperature in Figs. [11] and [T2j The diquark masses equal the pion mass at 
fi=0 due to the Pauli-Giirsey symmetry. At fi > one of the diquark mass 
decreases linearly and the other one increases linearly with \i as discussed 
in Appendix B. This is due to the trivial fact that there exist diquark and 
antidiquark states with energies at rest Ed = ±M<j = ±140 MeV. Their 
excitation energies change linearly with = 2/i, but with opposite signs. 

After diquark condensation, which occurs at /i = -m^ (fi B = m n ), the 
pion mass increases linearly with [i as discussed in Appendix B. As for the 
sigma and diquark masses, they couple after diquark condensation, |A | ^ 0. 
The mass matrix is 

fq 2 -m 2 -2A*m -2 Am \ 

M(uj,q)= -2Am \(q 2 - m 2 d ) + 2/jw -A 2 . (71) 

\-2A*m -A* 2 \{q 2 - m 2 ) - 2/iu ) 

where we have written q 2 = u> 2 — q 2 . The mass becomes 

m 2 = (M 2 - 2s 2 / 2 + 2s 2 / |A| 2 )/G ?0 2 / ) = 2|A| 2 . (72) 

The mass spectra for sigma, diquark and antidiquark are obtained by solving 
the dispersion relation detM(o;,0) = with respect to u. One solution is 
bj 2 = 0, and we are left with a 2nd order equation for u 2 , which is solved 
analytically as 

u 2 = 2|A| 2 + 10/i 2 + 2m 2 ± y/(6fj? + 2|A| 2 + 2m 2 ) 2 - 48/i 2 m 2 . (73) 

We show in Fig. [12] the results of the diagonalized masses. The coupling of 
the sigma meson and the diquark-baryons is very large. When the coupling 
is suppressed the sigma mass drops slightly as /i approaches the crossover 
chemical potential /i ~ 150 MeV and then increases together with the pion 
mass as shown in Fig. [TT] At the same time, the antidiquark mass in- 
creases rapidly with the chemical potential after diquark condensation. On 
the other hand, in the case of strong scalar-diquark coupling (Fig. [T2|) . the 
quantum-mechanical "non-crossing" rule is at work. The sigma mass in- 
creases continuously and the antidiquark mass joins the slowly increasing 
pion mass. 



24 



4-5. Equation of state of quark-hadron matter 

Next, consider the equation of state (EOS) of quark-hadron matter at 
various chemical potentials [i as functions of the temperature T. The pres- 
sure is given by the thermodynamical potential, p = —Q. For the case of 
only quarks, Q = &mf ( fI7|) . the pressure is essentially zero as shown by the 
thin solid curve in the confined region (small temperature), and gradually 
increases with temperature as deconfinement sets in featuring the entangle- 
ment of chiral symmetry and Polyakov loop effects. Finally the pressure 
becomes large at high temperature above the crossover transition. 
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Figure 13: The pressure divided by 
the Stefan-Boltzmann pressure for 
various chemical potentials \i =0 
(solid), 30 (dashed), 60 (dash-dotted) 
and 66 MeV (dotted) below diquark 
condensation as functions of temper- 
ature. Shown also is the one with 
only the quark degree of freedom 
(thin solid curve). 
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Figure 14: The logarithm of the pres- 
sure log(P) for various chemical po- 
tentials for [i = Q, 60 MeV below 
H c = 70 MeV and for [i = 80, 100 
and 200 MeV above [i c as functions 
of temperature. The unit of pressure 
is GeV 4 . 



In the confined region, the essential active degrees of freedom are hadrons. 
Their dynamics is governed by the hadron Lagrangian derived previously, in 
which mesons and diquark-baryons interact. We calculate thermodynamical 
potential, taking only the mass terms and integrating out the hadron fields 
by the Gaussian approximation as explained in Sect. 3. As mentioned pre- 
viously, the divergent zero-point energy terms in Eq. ( 165]) are regularized 
by introducing the cut-off momentum A. The zero point energy at T = 
and [i = is then subtracted from the thermodynamical potential. For the 
case A = 0, this procedure makes the contribution of the zero-point energy 
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terms vanish. Hence, the whole contribution to the pressure from the hadrons 
comes from the temperature dependent terms for A = 0. 

Shown in Fig. [13] is the pressure at fi = 0, 30, 60 and 66MeV. Given the 
small pion and diquark-baryon masses, the contributions of these degrees 
of freedom make the pressure significantly different from that with quark 
degrees of freedom only. This result is qualitatively similar to the one of 
the color SU(3) PNJL model |l4j, but here, due to the additional pressure 
of the diquark-baryon fields, the effect of the hadron contributions is much 
larger. The pressure increases rapidly as the chemical potential approaches 
the critical chemical potential fi c = This rapid increase of the pressure 

is caused by the diquark-baryon mode whose energy drops as ujd = — 2ji 
with increasing /i. When the temperature becomes small (T ~ lOMeV), we 
can expand the logarithm and the pressure can be written approximately as 

p~T [ _^L e -K®-2 M )/T _ T f qmax 5^9 e -K(S)-2,i)/T ^ ( 74 ) 



(27r)3 J q 27T 

where q max is appreciable only for the diquark-baryon mode (would be zero 
mode) as /i — )■ \i c . Since the pressure is divided by the Stefan-Boltzmann 
pressure, the pressure ratio close to the critical temperature shows rapid 
growth when /i approaches the critical chemical potential. 

Consider now the pressure over a wide range of chemical potentials /j,, 
in particular above the critical fi c for diquark condensation. The pressure is 
then dominated by quarks through their zero point motion renormalized by 
the vacuum value, with additional contributions from the zero point motion 
of mesons and diquark-baryons. Given that these zero point motion effects 
do not vanish at zero temperature, it is more appropriate to present the pres- 
sure as such, not divided by the Stefan-Boltzmann pressure, on a logarithmic 
scale. Results are shown in Fig. HH The zero point motion is largely influ- 
enced by the presence of the Bogoliubov spectrum of the diquark zero mode, 
discussed in Appendix B jgl, 0, 41, 48]. The pressure curves in Fig. [T41 are 



displayed in the whole fi range as functions of temperature. Shown in Fig. 
[H] are the pressures at /i =0, 60 MeV below /i c = m 7r /2 and \i = 80, 100 and 
200 MeV above fi c . The pressure below the critical chemical potential drops 
to zero at T = 0, while the pressure at \i > \i c stays finite at zero temper- 
ature. The pressure at low temperature increases rapidly across the critical 
chemical potential. The pressure at high temperature is dominated by the 
de-confined quark contribution and insensitive to the chemical potential. 
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4-6. Baryon density 

An interesting point of comparison with color SU(2) lattice simulations 
concerns the quark density as a function of chemical potential /i at zero tem- 
perature. When the baryon number symmetry is broken, the diquark con- 
densate becomes finite and the diquark-baryon becomes a Nambu-Goldstone 
boson. From the onset of diquark condensation, the quark number becomes 
finite. The quark density derived from the thermodynamical potential at 
mean field level is: 

Here, we write only the zero point oscillation terms, dropping the tempera- 
ture dependent terms. 




Figure 15: The quark density in unit of fm~ 3 as a function of the chemical potential 
[i in MeV. Shown by various points are the results of SU(2) lattice simulations [lil ]. 
The dashed curve denotes the result of the mean field approximation and the solid 
curve the results of the mean field and Gaussian approximation of the hadron 
contributions. 



In Fig. [15] the quark density is presented as a function of the chemical 
potential /z. As compared with lattice simulations [l9[ the mean field quark 
density comes out to be smaller than the lattice results by about a factor 
of two. The quark density becomes non-zero when the diquark condensate 
develops, starting at /i = \ra 1 ,. The quark density increases with the chemical 
potential. 
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Going beyond mean field level implies adding hadronic contributions to 
the quark density. In a first step we use the Gaussian approximation, drop- 
ping higher order terms of the hadron Lagrangian, and integrate out the 
meson and diquark-baryon fields. The hadronic part of the density is found 
by taking the derivative of flhadron, Eq. (JB"5]) . with respect to the chemical 
potential \x: 



Phadrc 



hadron 







dp, 



d 3 q 

dfi J (2^f 



^uJi(p) 3 



(76) 



where Ui (i = a, d, d*) and u„ are again the scalar boson, diquark, antidi- 
quark and pion energies, respectively. Only the temperature independent 
terms are written. Note that phadron vanishes for A = 0, while it becomes 
finite in the diquark condensed phase. The zero point energies are calcu- 
lated numerically and the momentum integrals are performed introducing 
the NJL cut-off A = 657 MeV. We varied the cut-off momentum to 2 GeV, 
and checked that the density increases only slightly. The result is shown in 
Fig. [15] by the solid curve. Evidently, using the Gaussian approximation, 
the effect of the hadron fields on the quark density is very small. This was 
anticipated by He 37j in their analysis of the hadron contributions. 

As already mentioned, the Gaussian approximation misses important 
hadronic interaction terms generated by the higher order pieces, and 
of the hadron effective Lagrangian. These interactions include, for ex- 
ample, scalar boson exchange between diquarks and mesons in various pos- 
sible combinations. The strength of these couplings is controled by constant 
A in Eq. ( I54p . We find A = 33 with the present parameter set. Altogether, 
the net attraction provided by such mechanisms is expected to increase the 
density significantly as a function of \x. A systematic treatment of these ef- 
fects is relegated to separate forthcoming work, using the Gaussian functional 



method discussed by Nakamura and Dmitrasinovic [49 



5. Summary and Conclusions 

In the present work we have investigated two-color quark-hadron matter 
at finite temperature and baryon chemical potential with the primary aim of 
studying, in a systematic and consistent way, the emergence and dynamics 
of baryonic degrees of freedom in addition to mesonic modes. Baryons are 
realized as diquarks in N c = 2 QCD. Their role as bosons is, of course, 
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fundamentally different from the N c = 3 case in which baryons are fermions 
subject to the Pauli exclusion principle. Nonetheless, the N c — 2 case is of 
conceptual interest because corresponding lattice simulations are not limited 
by the sign problem at non-zero, real chemical potential /x. 

The modeling is performed using the NJL approach supplemented by 
Polyakov loop dynamics. A full bosonization leads to an extended linear 
sigma model incorporating diquark degrees of freedom in the effective La- 
grangian. The thermodynamical potential derived from this effective La- 
grangian is then studied in the mean field approximation and beyond. Hadron 
properties are worked out as functions of temperature and chemical potential, 
and the phase diagram of quark-hadron matter is examined. The behavior of 
chiral and diquark condensates and their intertwining at finite /x, reflecting 
the underlying Pauli-Giirsey symmetry, is discussed. 

The Polyakov loop plays an important role at finite temperature. The 
deconfinement transition is now correlated with the chiral condensate. The 
characteristic temperatures for the chiral (T x ) and deconfinement (T dec ) cross- 
overs become about equal ( T x ~ T dec ~ 225 MeV) once the quarks are cou- 
pled to the Polyakov loop. At finite chemical potential and low temperature, 
T < T^c, the diquark condensate plays an important role. A non-zero di- 
quark condensate emerges at /x — \m 7T , the onset of Bose-Einstein condensa- 
tion in this model. As /x increases the diquark condensate grows quickly and 
stays finite at large chemical potential. With increasing diquark condensate, 
the chiral condensate decreases subject to the condition <j\ + |A | 2 — const. 

The meson and diquark masses behave naturally in accordance with the 
symmetry breaking pattern associated with chiral and diquark condensation. 
The pion mass first stays constant and then grows linearly with the chemical 
potential once diquark condensation starts. The sigma boson mass is about 
twice the dynamical quark mass until the onset of diquark condensation. 
Due to the strong coupling of the sigma boson with diquark-baryons, the 
mixed sigma-diquark modes behave quite differently from the case without 
coupling. 

The equation-of-state (EOS) of the quark-hadron system has interesting 
properties as a function of increasing chemical potential. When /x approaches 
its critical value /x c = 1/2 m, from below, the pressure divided by its Stefan- 
Boltzmann limit increase rapidly at low temperature. This reflects the fact 
that the energy of the lowest diquark mode drops as m n — 2fi and hence 
contributes prominently to the pressure. It becomes a zero mode at /x = /x c , 
the onset of diquark condensation, at which point the system develops a Bose- 
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Einstein condensation (BEC) phase. As the chemical potential increases 
beyond // ~ m n , the dynamical mass of quark quasiparticles drops below fi 
and the system undergoes a BEC-BCS crossover. With p further increasing, 
the quark quasiparticle energy develops a minimum at finite momentum p. 

Above the deconfinement transition, the EOS is governed by quark de- 
grees of freedom, but the hadronic modes still have a significant influence up 
to the crossover temperature range. 

The quark density p as function of the chemical potential is a quantity 
of interest that is readily accessible in two-color lattice QCD. It turns out 
that the mean field approximation of the (P)NJL model underestimates this 
density by about a factor of two. Corrections treated in Gaussian approx- 
imation do not change this result significantly. This is not surprising since 
the Gaussian approximation misses important correlations between diquarks 
and scalar bosons generated by the higher order terms (£ < - 3 - ) and C^) of the 
hadronic effective Lagrangian. The net attractive interactions produced by 
strong couplings in and are expected to raise the density consider- 
ably. Work along these lines using Gaussian functional methods is in progress 
and will be reported elsewhere. Nonetheless, the present work establishes a 
baseline for the further systematic expansion of fluctuations and correlations 
beyond mean field. 
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Appendix A. Meson and diquark-baryon fields and propagators 

Generic terms emerging in the expansion of the logarithmic term tr ln(l + 
SK) = J2k are given in this Appendix. We write here the k = 2 term 
for the kinetic and mass terms, the k = 3 term for the coupling terms and 
the k = 4 term for the interaction terms. 
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Appendix A.l. Pion 

The result of the Fourier transformed propagator for the pion in Eq. (1381) 

is: 

T PP (g) = '-tr J -0^g 2 o (H; ll5 H+ l5 + H+ l5 H~ l5 (A.l) 

where p' denotes p + q with q the outer four momentum. We have to take 
care of the spinor trace and the divergent integrals I2 and Iq. For q = one 
finds: 

r„(0) = ^ — W - {E ^ M - Im (A.2) 



2 



(Po-E p -)( P o + E;)-2A 
(^-(^a) 2 )(^-(^) 2 ) 

= glh - 2gfa 2 io ■ 

The first order derivative term is zero and we work out the second order 
derivative term. After a lengthy manipulation we get 

d^Vppiq)]^ = g%I g^ ■ (A.3) 

Appendix A.2. Sigma boson 

Next we consider the a boson term: 

T ss (q) = ~tr J -J^o 2 (G+G+ + G p/ G p + H p ,H+ + H}H;) . (A.4) 

In the g = limit: 

T ss (0) = glh - 2^ 2 I - 2^ 2 m 2 / . (A.5) 

The first order derivative term vanishes. The result of the second order 
derivative of Y ss is 

d^d v T ss (q)\ q=0 = glhg^ u . (A. 6) 

The pion field does not mix with other fields due to its pseudoscalar nature. 
However, the sigma boson field mixes with diquark and antidiquark fields. 
The mixing terms will be calculated later. 
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Appendix A. 3. Diquark-baryon 
We calculate the dd* term 



F d d*(q) = |tr J j^g 2 d (-G^ 5 t2T 2 Gp^t 2 T2) . 
For q = 0, we get 

d 4 P ( Vl ~ {E-f pg - (25+) 



rv(o) = -^tr 



+ 



= \g 2 d h-\g d \^h. 

The first order time derivative is 



(2tt)4 V(p§-(^a) 2 ) 2 (Pl-(EI) 2 ) 2 

2 



(A.7) 



(A.8) 



<%r Q 



+ 



d A p 



E n 



(Pi - (EX) 2 ) 2 (P 2 - (£ A ) 2 ) 2 



(A.9) 



+ 



(p 2 o-(Et) 2 ) 2 (P 2 -(E A 



i2\2 



= 9&Io ■ 
The second order derivative of Thh* is 



1 



Next, the d*d term is: 



r d . d (g) = -#jtr 



(2^ 



'-G , p ,75^2r 2 G'+75t2r2) . 



For g = 0, we find 



r^(0) = -gjl 2 - -gl\Ao\% . 

The first order time derivative becomes: 

d T d * d (q)\ q=0 = -gjfilo . 
The second order derivative term is: 

d„d u r d * d (q)\ q=0 = ^gjlog^u ■ 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



(A.14) 
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Now we consider the dd and d*d* terms that describe diquark mixing. The 
dd term is: 

j^g 2 d H+ l5 t 2 r 2 H+ l5 t 2 r 2 . (A.15) 

For q = 0, we get 



U, y (27r)4 V(p 2_ (E + )2)2 (pg _ (£^2)2 

= -^S(AS) 2 / . (A.16) 
All the derivative terms are zero. By analogy the d*d* term becomes: 

IW(g) = |tr y j^g 2 d H-,-f 5 t 2 T 2 H~-f 5 t 2 T 2 , (A.17) 
iW(o) = -^(A ) 2 / . 

Appendix A. 4- Sigma- diquark mixing 

The pion does not mix with diquark-baryon fields due to its pseudo-scalar 
nature. However, the sigma boson mixes with diquark fields when diquark 
condensation takes place. We calculate the mixing terms sd, ds, sd* and d*s. 

T ds (q) = l -trj ^(-g g d )( G ; /l5 t 2 r 2 H+ + H} l5 t 2 r 2 G;) .(A.18) 

For q = 0, we get 

r ds (0) = -(? (? 2 A*m/o . (A.19) 

For the sd term, 

^(-^)(iJ+G p + 75t 2 r 2 + G pl H; l5 t 2 r 2 ) . (A.20) 

At g = 0: 

r sd (0) = -^o^A*m7 • (A.21) 
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The derivatives of T^s and T s d terms are finite. We calculate the d*s term 

j^- 4 gMG pll5 t 2 T 2 H- + H;, l5 t 2 r 2 G;) . (A.22) 

For q — 0, the result is: 

rV s (0) = -g g 2 d A mI . (A.23) 

Next, the sd* term: 

/d^ 

j^9o9d{H pl G pl5 t 2 r 2 + G+H pl5 t 2 r 2 ) . (A.24) 

Its q = limit is: 

r sd *(0) = -g g 2 d A mI . (A.25) 
All k = 2 terms are collected in Eq. (|4T)|) . 

Appendix A. 5. The third order term 

Here we examine = —^tr(SK). This term generates leading com- 
binations of G ± and H ± with 75. With the principle of taking only the 
divergent integrals, the possible combinations are GGG or GGH with 75. 
From dimensional analysis, the coupling constant should have mass dimen- 
sion one. Thus for baryon-number conserving channels, the relevant factor 
is m. For baryon-number non-conserving channels the factors are A or 
Aq. Hence, the terms for baryon-number conserving channels are generated 
from GGG and the terms for baryon-number non- conserving channels are 
generated from GGH. We find 32 terms that remain finite from the above 
discussion. We do not write all the terms, but give a representative example 
for a baryon-number conserving channel: 

U$> = ^ 3 tr/-|^(G + G + G + + G-G-G-) S 3 (A.26) 
= -2^/ ms 3 , 
and for a baryon-number non-conserving channel: 

U% = ^tr|^(G + 75^' + 75G-75 + ^' + 75G-7 5 G + 75(A.27) 
+G- l5 G + l5 H' + l5 )d 2 d* = -g 3 d I A* d 2 d* , 
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where 

**~i£ggr>+-ir&r h - (A - 28) 

The final result for k = 3 is displayed in Eq. f l42|) . 

Appendix A. 6. The fourth order term U^' 

Next, consider = ^tr(SK) 4 . By examination it turns out that the 
only possible combination is GGGG and the coupling constants are dimen- 
sionless. Out of the many terms that appear we give here a generic example: 

UP = ^o 4 tr|^(G + 75G + 75G + 75G + 7 5 + (+^-)K 4 .(A.29) 

The products of gamma matrices are worked out as follows: 

A_7o75A+7o75A_7o75A+7o75 = A_ , (A.30) 
A+7o75A_7o75A+7o75A_7o75 = A + , 

and all other terms vanish. From the first term one finds: 

tr , G+75 « +75G+75G+75 = 2 (j^^ + ^ZJ^y) ■ (A.3D 

The factor 2 in this equation comes from the spinor trace. The same result 
is obtained from the second term of Eq. ( 1A.29I) . hence the result of the 
complete 7? 4 term is 

i 4 ^ f d^p 



9>f* I 77^ ( -2 rE^ + -2 rE^vt ) < A-32 ! 
1 



2 yu r °J (2tt) 4 V(^-(^a) 2 ) 2 (pg-(JSS wa 



- 4 t -*4 

Consider another example, the combination involving 7f 2 dd*, which has 12 
terms of the form GGGG. By repeating the above procedure, one finds: 

{ (A-kn + (ri-(W)O ^ (A ' 33) 

= -ghlh^dd* . 

The complete result for = 4 is written in Eq. ( jUj) . 
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Appendix B. Hadron masses 

Here we derive the mass spectrum of hadrons. The pion does not mix 
with other hadrons, but the sigma boson mixes with diquark-baryons once 
diquark condensation takes place. We list all the masses of mesons and 
diquark baryons. 

Appendix B.l. Pion and sigma masses 

The renormalized pion mass is written as 

ml = (M 2 - 2g 2 I 2 + Ag 2 I^ 2 )/(g 2 I ) . (B.l) 

From hereon, all the masses are understood to be renormalized ones. Con- 
sider the case of |A | = and define divergent integrals for // = as: 

J 2 for finite /i with |A | = is written in terms of I 2 and Jq 

h = rtr * J (pTW? + W^f) (B ' 3) 

= /? + 2«. 
In the region of interest, I = Iq and m 2 becomes 

ml = (M s 2 - 2glll)/{glll) = m 2 . (B.4) 

Hence, the pion mass does not depend on /i for |A | = 0. 

The diquark gap equation with the Pauli-Giirsey symmetry G = H 
provides g = g d and M s = M d and therefore M 2 S — 2g%I 2 = for the finite 
|A | in which ml = 4fi 2 . This means that the pion mass in the diquark 
condensed phase increases linearly with //. 

As for the sigma mass, m 2 s is larger than ml by twice dynamical quark 
mass squared, 4m 2 , for the case without sigma-diquark mixing: 

m 2 s = (M 2 - 2g 2 I 2 + 4g 2 I of i 2 + 4g 2 I m 2 ) / (g 2 I ) = m 2 + 4m 2 . (B.5) 

Hence as the dynamical quark mass approaches zero, these two masses coin- 
cide. The sigma boson, however, mixes with diquark-baryons once diquark 
condensation sets in. 
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Appendix B.2. Diagonalization of sigma and diquark-bary on mass matrix 

The sigma boson mixes with diquark-baryons, and furthermore diquarks 
mix with each other. The sigma and diquark Lagrangian density with kinetic 
and mass terms is written in matrix form: 

1 fd 2 +m 2 s 2A*m 2 Am \ 
Cf d = --$M 2Am \{d 2 + m 2 d ) + 2ifi8P A 2 $ , 

2 \ 2A*m A* 2 \{d 2 + m 2 ) - 2ifxd°J 

(B.6) 

with the scalar fields representation $* = (s, d, d*). The diquark-baryon mass 
becomes: 

m 2 d = (M 2 -2g 2 I 2 + 2g 2 I \A\ 2 )/(g 2 I ) (B.7) 
m 2 -4/i 2 |A| = , 
2|A| 2 |A|^0, 

using the diquark gap equation Mj — 2(yf^/ 2 = and ga = g for |A| ^ 0. 

We write the mass matrix in the momentum representation (Minkowski 
space) as 

V - m 2 -2A*m -2Am \ 

M(u, q) = ( -2Am i(g 2 - m 2 ) + -A 2 , (B.8) 

-2A*m -A* 2 \(q 2 - m 2 ) - / 

where q 2 = u 2 — q 2 . The mass spectra for sigma, diquark and antidiquark 
are obtained solving the dispersion relation det M(u, 0) = with respect to 
bj 2 . One of the solutions is zero and identified with the Nambu-Goldstone 
(NG) boson at |A| ^ 0, since the baryon number symmetry is spontaneously 
broken. The other solutions are given analytically as 



u 2 



2|A| 2 + 10/j, 2 + 2m 2 ± v/(6/x 2 + 2|A| 2 + 2m 2 ) 2 - . (B.9) 



In the limit of |A| — > at the onset of diquark condensation (/x = m^/2) the 
dispersion Eq. (1B.9|) becomes 
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The upper solution is the squared sigma mass m 2 (without mixing effect) 
and the lower one is the squared antidiquark mass at /i = m n /2. In the limit 
m — > at large \x the dispersion Eq. (IB.9j) gives 



- 2 = V 2 =m '- 2 (B.H) 
* 16/x 2 + 4|A| 2 



One solution represents the linearly increasing pion mass. The other one is 
the squared antidiquark baryon mass at large /i where mixing effects with 
the scalar boson disappears. 

When A = (/i < m n /2) the mixing terms are zero. The diquark masses 
are obtained as 

u = m n ±2fi. (B.I2) 

The diquark and antidiquark masses move linearly with for A = 0. The 
diquark-baryon behaves as m n — 2[i and becomes the NG boson at non-zero 
|A|. On the other hand, the antidiquark with baryon number —1 has the 
mass m n + 2/i. These mass spectra including the pion are plotted in Fig. [I2J 
In the absence of the sigma-diquark mixing term, — 2mA ( -*\ the diquark 
masses become 

u A . = v/l6/i 2 + 4|A| 2 , (B.13) 
oo A = . 

The sigma mass approaches the pion mass as the dynamical quark mass tends 
to zero due to the chiral symmetry restoration. This feature is shown in Fig. 

o 

Appendix B.3. The Bogoliubov mode 

The hadronic contributions to the pressure and the quark density in- 
volve spectra of hadrons. It is important to know how the zero mode be- 
haves at finite momentum in the diquark-condensed phase. The Bogoliubov 
excitation jsl, 0, 41] for the zero mode results from the solution u(q) of 



det M(u, q) = in Eq. flB.8j) at finite momentum q as 



det M{u,q) = tu 6 - [{q 2 + m 2 ) + 2(q 2 + m\ + 8/i 2 )]w 4 (B.14) 
+ [q 2 (q 2 + 2m 2 ,) + 2(g 2 + m 2 )(g 2 + m 2 d + 8fi 2 ) - 16m 2 |A| 2 ]w 2 
-q 2 [{q 2 + m 2 ){q 2 + 2m 2 ) + 16m 2 |A| 2 ] = . 
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We write the solutions as Ui (i = s, d, d*). Dropping the sigma-diquark 
mixing terms in the full mass matrix ( 1B.8|) . one arrives at an analytical 
expression for oo(q) of the Bogoliubov mode. Consider the reduced mass 
matrix 

n , / u 2 - q 2 - m 2 + AfjM -2A 2 \ 

D ( u > $ = { -2A* 2 co 2 - e -ml- 4^J ' (B - 15) 

and solve the dispersion relation det D(u,q) = with respect to w(q) for 
the lowest mode: 

u 2 = g 2 + 2|A| 2 + 8/i 2 (B.16) 

-V(q 2 + 2|A| 2 + 8/i 2 ) 2 - q 2 (q 2 + 4|A| 2 ) 
g 2 (g 2 + 4|A| 2 ) 



2(g 2 + 2|A| 2 + 8/i 2 



In the BEC limit, the pion mass m n = 2/x is much larger than |A|. The 
Bogoliubov excitation is then written as 




^>1^he + — • (B.17) 



The zero mode varies linearly with \q\ at small momentum. With inclusion 
of the sigma-diquark coupling, the dispersion equation for oo(q) with the full 
mass matrix is solved numerically. For m n 3> |A| one confirms that the zero 
mode has the Bogoliubov excitation spectrum as expressed in Eq. (1B.17j) . 
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